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ON  THE  SIMULATION  OF  TRANSONIC  SHOCK-TURBULENT 
BOUNDARY  LAYER  INTERACTIONS  IN  CRYOGENIC  OR 
HEAVY  GAS  WIND  TUNNELSt 


G.  R.  Inger* 

Virginia  Polytechnic  Institute  and  State  University 
Blacksburg,  VA 


Abstract 

The  role  of  the  basic  similitude  parameters  governing  transonic  normal 
shock- turbulent  boundary  layer  interaction  effects  in  cryogenic  wind  tunnel 
tests  is  studied  theoretically  for  the  non-separating  case.  Besides  Mach  and 
Reynolds  number,  these  parameters  are  the  wall  to  total  temperature  ratio, 
specific  heat  ratio  y»  viscosity- temperature  exponent  and  Prandtl  number.  The 
results  show  that  lack  of  temperature  ratio  simulation  has  a significantly 
adverse  effect  on  interactive  skin  friction  and  hence  separation  onset  compared 
to  the  adiabatic  free  flight  case;  higher  y's  than  air  also  may  have  some 
effect. 


t Based  on  work  supported  by  the  Office  of  Naval  Research  under  Contract 
N00014-75-C-0456. 

* Professor  of  Aerospace  and  Ocean  Engineering,  Assoc.  Fellow,  AIAA 
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NOMENCLATURE 

skin  friction  coefficient 

pressure  coefficient  [2(p-Po)/pQgUQg2] 

boundary  layer  profile  shape  factor  (6*/e*) 

= = Mach  number  of  incoming  in viscid  flow 

static  pressure 
Prandtl  number 

Reynolds  number  based  on  and  L 

recovery  factor 

absolute  static  temperature 

velocity  components  in  x and  y directions,  respectively 
coordinates  parallel  and  normal  to  surface 


Y ratio  of  specific  heats 

6 boundary  layer  thickness 

6*  boundary  layer  displacement  thickness 

aP  static  pressure  jump  across  incident  shock 

n non-dimensional  displacement  surface  perturbation  (Ay/s^) 

y coefficient  of  viscosity 

u)  viscosity- temperature  exponent 

p density 

e*  boundary  layer  momentum  thickness 

Superscri pts 

( )'  disturbance  quantity  (with  respect  to  incoming  flow) 
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Subscripts 

e freestream  properties  upsteam  of  shock  at  boundary  layer  edge 

0 denotes  undisturbed  (not  stagnation)  flow  property  ahead  of  shock 

L length  of  undisturbed  shock  location 

w property  at  wall  value 

i - 1,  2,  3,  denotes  region  of  variable  (Fig.  1) 
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INTRODUCTION 


I 


It  is  well  known  that  the  aerodynamics  of  missiles,  entry  vehicles  and 
aircraft  flying  at  transonic  speeds  is  significantly  influenced  by  the  inter- 
action between  local  normal  shock  waves  and  the  viscous  boundary  layer  flow. 
Consequently,  it  is  important  to  insure  that  this  flow  feature  is  adequately 
simulated  in  experiments  carried  out  in  future  cryogenic  or  heavy  gas  high 
Reynolds  number  wind  tunnels,  especially  since  their  cost  per  data  point  is 

expected  to  be  higher  than  usual. ^ In  addition  to  the  flight  Mach  and  Reynolds 

1 2 

numbers  which  are  by  design  simulated  in  such  facilities,  ’ four  other  para- 
meters are  involved  in  the  interaction  effect  on  a transonic  flow  field  which 
may  not  be  duplicated  owing  to  the  very  low  temperature  -high  pressure  working 

fluid  involved:  wall  to  total  temperature  ratio  T /T. , specific  heat  ratio  y, 

w t 

viscosity  temperature  dependence  exponent  a)(p-T“)  and  the  Prandtl  number  Pr. 

Of  these,  the  first  is  deemed  especially  important  since  the  aerodynamic  tests 
in  some  proposed  short  duration  cryogenic  transonic  wind  tunnels  result  in  the 
model  being  at  much  higher  temperature  than  the  flow  total  temperature  during 
much  of  the  test.  Moreover,  the  specific  heat  ratio  of  cryogenic  nitrogen  flow 
can  be  significantly  higher  (y-1.5-1.8)  than  air,  as  can  that  of  argon 
(y  = 1.67),  and  this  may  have  some  influence  on  the  interaction  effects  since 
Y is  involved  in  both  the  shock  and  boundary  layer  governing  equations. 

Specific  heat  ratios  lower  than  air  (y^I.1)  are  also  of  interest  in  facili- 

4 

ties  using  heavy  gases  such  as  Freon  12. 

Although  the  influence  of  real  gas  effects  on  the  inviscid  transonic 
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flow  ’ and  on  laminar  boundary  layer- supersonic  oblique  shock  interaction 

has  been  examined  in  detail,  an  appraisal  of  the  more  important  turbulent- 


normal  shock  interaction  problem  has  yet  to  be  made.  Existing  theoretical 


j studies  of  transonic  shock  turbulent  boundary  layer  interaction  to  date  have 

1 

been  confined  to  the  case  of  adiabatic  flow  in  normal  air.  To  appraise  the 
influence  of  the  aforementioned  similitude  parameters,  the  present  paper  there- 
fore describes  the  extension  and  application  of  a recently- developed  basic 
analytical  theory  of  the  nonseparating  transonic  interaction  of  a weak  normal 
shock  with  a turbulent  boundary  layer  to  include  the  effect  of  heat  transfer 
and  arbitrary  values  of  y,  oi  and  Pr.  Results  are  given  for  pressure  distri- 
butions, flow  geometry  and  skin  friction. 

OUTLINE  OF  THE  THEORETICAL  MODEL 

For  unseparated  turbulent  interacting  boundary  layers  (M^  5 1*3)  it  is 
possible  to  construct  a fundamentally-based  theory  of  the  problem.  In  its 
original  form  the  flow  is  taken  to  consist  of  a known  adiabatic  zero  pressure 
gradient  boundary  layer  profile  M^Cy)  subject  to  small  transonic  disturbances 
due  to  an  impinging  weak  normal  shock.  The  theoretical  model  of  this  inter- 
action emerges  from  an  asymptotic  analysis  of  the  compressible  Navier-Stokes 
equations  at  High  Reynolds  numbers,  giving  a linearized  inviscid  boundary 
value  problem  surrounding  the  nonlinear  shock  discontinuity  and  underlaid  by  a 
thin  viscous  disturbance  sublayer  as  schematically  illustrated  in  Figure  1.  As 
described  in  detail  elsewhere,®  an  approximate  analytical  solution  can  be 
achieved  in  the  leading  approximation  by  introducing  suitable  simplifying 
assumptions  which  are  physically  sound  provided  = M^^  is  not  too  close  to 
unity  (M.|  js  1.05);  solutions  can  then  be  obtained  by  operational  methods  for 
each  subregion  indicated  in  Fig.  1.  An  arbitrary  value  of  y may  be  retained 
in  so  doing.  Detailed  comparisons®*^  have  shown  that  the  results  give  a good 
account  of  the  important  engineering  aspects  of  the  interaction  (e.g.,  ap/ay 
effects  across  the  boundary  layer  near  the  shock,  interaction  pressure  rise 
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along  the  wall,  displacement  thickness  growth  and  interactive  skin  friction) 
over  a wide  range  of  Mach- Reynolds  number  conditions  at  very  low  computational 
cost.  It  is  also  readily  adapted  to  non-adiabatic  flow  conditions  with 
arbitrary  values  of  to  and  Pr,  and  hence  is  well -suited  as  a basis  for  the 
present  study. 

The  matching  of  the  aforementioned  regional  solutions  yields  readily- 

solved  linear  integral  equations  for  the  disturbance  pressure  along  both  the 

boundary  layer  edge  and  wall.®  The  remaining  interactive  flow  properties  can 

then  be  determined  in  terms  of  p'.  For  example,  double  integration  of  the 

w 

combined  continuity  and  y-momentum  equations  across  the  width  of  the  inviscid 
rotational  disturbance  region  in  the  boundary  layer  yields  the  local  increase 
in  displacement  thickness.  The  viscous  disturbance- flow  solution  in  the  under- 
lying viscous  sublayer  determines  the  interactive  skin  friction  perturbation 
associated  with  p'.  Thus,  upon  correcting  for  the  important  non-linear  inertia 

W 

effects  in  an  adverse  pressure  gradient  using  the  general  non-dimensional  wall 

O 

shear-pressure  solution  ahead  of  separation  given  by  the  triple  deck  theory 


(converted  to  turbulent  flow  by  expressing  all  results  in  terms  of  instead 

0 


where  the  non-dimensional  function  F is  essentially  unity  ahead  of  the  shock 

X < 0 and  vanishes  slowly  behind  it  with  F (x/6q)~^^^  far  downstream.  It 

is  seen  from  Eq.  (1)  that  depending  on  Reynolds  number  (C^  ),  a sufficiently 

0 

strong  interactive  pressure  rise  can  cause  incipient  separation  (C^  -*■  0)  near 
the  shock. 

9 

It  can  be  shown  that  the  presence  of  small -to  moderate  heat  transfer 


in  transonic  flow  does  not  significantly  influence  either  the  inviscid  or 
viscous  disturbance  fields  by  introducing  any  new  perturbation  equation  terms; 
in  the  leading  approximation  it  enters  only  implicitly  as  it  influences  the 
undisturbed  (y),  skin  friction  and  boundary  layer  thickness  6^.  The 

boundary  layer  profile  used  is  based  on  an  accurate  flat  plate  turbulent  eddy 

viscosity  model  with  arbitrary  w that  is  well-suited  for  interaction  studies^^; 

it  uses  a modified  Crocco  integral  to  compute  the  temperature  and  hence 

provides  for  an  arbitrary  degree  of  surface  heat  transfer  with  a recovery 
1/3 

factor  r *c  ' and  arbitrary  y 

It  is  noted  that  the  purpose  of  this  study  is  to  delineate  the  overall 
engineering  parameter  sensitivies;  accordingly,  we  have  deliberately  used  a 
very  simplified  real  gas  model  involving  constant  equivalent  thermodynamic 
(y)  and  transport  (u,  Pr)  properties.  Should  their  effect  be  significant,  it 
is  understood  that  a more  detailed  treatment  of  the  thermodynamic  and  transport 
relations  within  the  interaction  flow  analysis  might  be  necessary. 

DISCUSSION  OF  RESULTS 

To  provide  a general  basis  for  appraising  the  similitude  requirements  of 
various  facilities  with  different  gases  and/or  thermal  histories  (ranging 
from  fan-driven  or  blow-down  to  short  duration  Ludwieg  tubes),  a wide  range  of 
the  parameters  was  studied.  The  following  are  typical  results. 

Heat  Transfer  Effect 

The  predicted  wall  temperature  effect  on  the  interaction  pressure 
distributions  along  the  wall  is  illustrated  in  Fig.  2 and  is  seen  to  be  quite 
weak;  as  expected,  this  was  found  to  be  true  over  a range  of  shock  strengths 
and  Reynolds  numbers.  Increasing  wall  temperature  tends  to  slightly  increase 

the  upstream  influence  and  lower  the  pressure  immediately  downstream  of  the 

« 
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shock  with  a tendency  to  lengthen  the  asymptotic  "tail".  Regarding  the 
important  property  of  the  characteristic  upstream  influence  distance  x 

up 

(here  defined  as  the  distance  upstream  of  the  shock  where  the  local  inter- 

Q 

action-induced  pressure  rise  is  only  5%  of  the  overall  total),  we  have  found 
that  over  a wide  range  of  conditions  the  ratio  «.0(1)  and  is  only 

slightly  affected  by  heat  transfer;  thus,  e.g.,  cooling  reduces  the  upstream 
influence  essentially  proportional  to  the  corresponding  reduction  in  6^,  as 
also  observed  at  higher  supersonic  Mach  numbers. 

The  interaction-induced  growth  of  the  boundary  layer  displacement  thick- 
ness is  of  great  practical  interest  since  this  often  has  a significant  back- 
effect  on  the  in viscid  flow  and  shock  position  on  an  airfoil  or  in  channel 
flows. ^ Typical  results  are  shown  in  Fig.  3 illustrating  the  expected 
thinning  out  with  increasing  Reynolds  number  or  wall  cooling.  The  influence 
of  a hot  wall  T^  > T^  is  increasingly  significant  at  lower  Rej^.  Note  also 
that  this  effect  is  significant  in  terms  of  the  ratio  since  6^  is 

itself  increased  by  heating,  it  is  even  larger  in  terms  of  aY. 

The  effect  of  shock-boundary  interaction  on  the  local  skin  friction  is  of 
particular  importance  in  transonic  airfoil  design  and  testing,  since  it  bears 
directly  on  the  downstream  boundary  layer  behavior  and  its  possible  separation. 

Since  wall  temperature  influences  the  undisturbed  skin  friction  C,  (see 

'o 

Table  1),  the  relative  effect  on  its  interactive  decrement  alone  can  be  shown 

by  plotting  the  ratio  C^(x)/C^  in  Fig.  4A.  Owing  to  the  interaction-induced 

adverse  pressure  gradient,  C^/C^  typically  decreases  downstream  toward  the 

® * 

shock  with  a minimum  occurring  slightly  behind  it  , followed  by  a subsequent 

* Although  no  longer  valid  for  separated  flow  where  C^(x)  < 0 over  some  portion 
of  the  wall,  the  present  theory  is  still  useful  to  indicate  trends  toward  this 
situation,  i.e.,  where  and  when  incipient  separation  at  some  x)  first 

7 9 ■ 

occurs.  ’ 
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gradual  rise  further  downstream.  It  is  seen  that  wall  temperatures  appreciably 
different  from  adiabatic  significantly  influence  the  interactive-reduction  of 
and  hence  the  onset  of  incipient  separation;  in  particular,  hot  walls  are 
predicted  to  magnify  the  adverse  effect  of  the  interaction  of  and  hasten 
the  occurance  of  incipient  seoaration  under  the  shock,  whereas  wall  cooling 
has  the  opposite  beneficial  effect.  Judging  by  comparison  with  calculations 
showing  the  effect  of  Reynolds  number  (shown  in  Fig.  4B),  it  would  appear 
that  proper  wall  temperature  simulation  may  be  of  comparable  importance  to 
Reynolds  number  as  regards  skin  friction. 

The  foregoing  results  are  in  qualitative  agreement  with  experimental  data 
on  non-adiabatic  interactions  at  supersonic  speeds  with  oblique  shocks. 

However,  to  the  authors  knowledge,  there  exists  as  yet  no  experimental  data  on 
transonic  unseparated  interactions  with  heat  transfer. 

Specific  Heat  Ratio 

Although  studies  have  shown  that  the  purely  inviscid  aspects  of  the 
transonic  flow  field  are  not  significantly  affected  by  the  value  of  y in 
the  range  1.1  - 1.8  (Ref.  3-5),  coupling  with  viscous  effects  as  occurs  in 
shock-boundary  layer  interaction  may  change  this  conclusion  and  hence  warrants 
investigation. 

The  results  of  our  calculations  show  that  there  is  only  a barely-discernable 
effect  of  Y on  the  interactive  wall  pressure  distribution  over  a wide  range  of 
Mach  and  Reynolds  numbers:  as  shown  in  Fig.  5,  increasing  y slightly  expands 
the  nondimensional  streamwise  scale  of  the  interaction.  Likewise,  there  is 
only  a very  small  change  in  the  interactive-thickening  of  the  boundary  layer 
(Fig.  6).  These  predicted  insignificant  effects  are  in  agreement  with  experi- 
mental observations^  and  similar  to  theoretical  results  for  the  laminar  oblique 
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shock  interaction  case. 

The  influence  of  y on  the  interactive  skin  friction,  however,  is  more 

interesting:  Fig,  7 shows  the  typical  effect  on  the  ratio  C^(x)/C^  . Increas- 

0 

ing  Y reduces  not  only  the  undistrubed  C,  but  also  even  more  the  interactive 

^0 

decrement  so  as  to  moderately  hasten  the  onset  of  incipient  separation.  The 
explanation  for  this  evidently  lies  in  the  noticeable  increase  of  the 
adiabatic  wall  temperature  with  y (see  Table  1);  increasing  y thus  has  an 
adverse  effect  on  qualitatively  equivalent  to  that  of  a higher  wall  tempera- 
ture. Combined  with  the  foregoing  wall  temperature  effect,  this  lack  of  air 
property  simulation  in  a cryogenic  facility  with  significantly  higher  y may  be 
important  where  skin  friction  and  incipient  separation  effects  are  significant, 
especially  in  view  of  the  drastic  change  in  the  basic  interaction  flow  pattern 
when  separation  occurs  (Fig.  8). 

Viscosity  Exponent  and  Prandtl  Number 

Examination  of  available  thermophysical  property  data  for  cryogenic 
12 

nitrogen  indicated  that  over  the  working  temperature  range  of  aerodynamic 

testing  interest,  a power  law  viscosity- temperature  relationship  y^T“  is  a 

reasonable  approximation,  with  u)  ss  .75  and  decreasing  slightly  with  decreasing 

absolute  temperature.  To  assess  the  sensitivity  of  the  viscous  interaction 

to  w,  calculations  were  made  for  w = .5  (classical  kinetic  theory),  .75  and 

.80  i the  results  showed  a completely  negli gable  effect  on  all  the  interaction 

properties  of  physical  interest  including  the  skin  friction. 

12 

Likewise,  Prandtl  number  data  shows  little  change  from  the  usual  air 
value  (x.72)  over  a wide  temperature  range  down  to  nearly  the  liquifaction 
limit;  since  the  influence  of  Pr  is  felt  through  the  basic  turbulent  boundary 
layer  recovery  factor  (r  as  Pr  ' ),  negli gable  influence  of  this  parameter  on 
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the  interactive  similitude  would  be  expected.  Calculations  in  which  r was 
varied  from  .8  to  1.0  verified  this. 

Upstream  Boundary  Layer  History 

In  addition  to  their  direct  effect  on  the  interaction,  the  aforementioned 

similitude  parameters  also  may  have  a significant  indirect  effect  through 

the  influence  on  the  upstream  boundary  layer  history  as  reflected  in  the 

1 3 

incoming  boundary  layer  profile  shape.  For  example,  it  is  wellknown  that 

shock-boundary  layer  interaction  can  be  very  sensitive  to  the  laminar-turbulent 

transition  history  including  tripping  devices  and  surface  roughness  (Fig.  9), 

and  these  in  turn  may  be  influenced  by  wall  heating  and  specific  heat  ratio 

effects. The  role  of  the  pressure  history  ahead  of  the  stock  may  also  be 

very  significant:  Fig.  10  illustrates  this  with  some  results  obtained  by 
1 8 

Panaras  and  Inger  that  show  the  appreciable  change  in  normal  shock-turbulent 
boundary  layer  interaction  properties  associated  with  small  variations  of  the 
initial  shape  factor.  Thus,  if  the  upstream  pressure  gradient  effect  on  the 
boundary  layer  (especially  if  adverse)  is  influenced  by  T and  y . as  the 

W ^ 

shape  factor  values  shown  in  Table  1 suggests  it  could  be  , this  must  be 
included  in  the  similude  appraisal  of  the  subsequent  interaction  region. 

CONCLUDING  REMARKS 

The  present  study  has  used  a basic  theoretical  model  as  a tool  to  appraise 
the  influence  of  nonadiabatic  wall  temperatures  and  lack  of  specific  heat  ratio 
simulation  on  the  shock-boundary  layer  interaction  aspect  of  the  transonic 
flight  regime  in  either  flight  or  wind  tunnel  applications.  For  example,  the 
sensitivity  of  skin  friction  predictions  to  T^/T.  and  y suggests  that  the 
lack  of  freeflight  adiabatic  wall  temperature  ratio  and  gas  property  simulation 
in  a cryogenic  tunnel  may  significantly  exaggerate  the  interaction  effect  and 
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its  attendant  flow  separation  on  the  model  compared  to  a flight  case  at  the 
same  M-j  and  Rej^.  Consequently,  some  basic  experimental  studies  of  transcnic 
interactions  under  non-adiabatic  conditions  appear  desireable  to  test  these 
conclusions. 

It  is  noteworthy  that,  aside  from  broadening  our  basic  understanding  of 
viscous-inviscid  interations,  there  are  two  other  practical  applications  where 
the  results  of  the  present  study  may  be  of  interest:  (a)  post-entry  transonic 
flight  phase  of  the  Space  Shuttle  orbiter,  where  transonic  shock-boundary 
layer  interactions  can  take  place  on  the  hot  surface  caused  by  the  entry 
heating  history;  (b)  studies  of  transonic  flow  around  cooled  turbine  blades 
operating  in  hot  gas  flow  conditions. 
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Temperature  Effect  on  Interactive  Thickening 
of  the  Boundary  Layer 


Temperature  Effect  B.  Reynolds  Number  Effect 


Fig.  5 


Specific  Heat  Ratio  Influence  on  Wall  Pressure 


Specific  Heat  Ratio  Effect  on  Interactive  Specific  Heat  Ratio  Effect  on  Interactive 

Thickening  Skin  Friction  Distributions 


pstream  Laminar-Turbulent  Transition  Effect 
on  the  Shock-Boundary  Layer  Interaction 


IKCUWTY  CLAaSIFICATlOW  OF  THIS  PAGE  (Whan  Dtt^Bnlmnd)^ 


REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

REPORT  NUMBER  2.  GOVT  ACCESSION  NO. 

VPI-AERO-080'^ 

S.  RECIPIENT’S  CATALOG  NUMBER 

4.  TITLE  fand  SuMtIm) 

ON  THE  SIMULATION  OF  TRANSONIC  SHOCK-TURBULENT 
BOUNDARY  LAYER  INTERACTIONS  IN  CRYOGENIC  OR 

HEAVY  GAS  WIND  TUNNELS 

S.  TYPE  OF  REPORT  S PERIOD  CDVERED 

S.  PERFDRMING  ORG.  REPORT  NUMBER 

7.  AUTHORfa; 

G.  R.  Inger 

6.  CONTRACT  OR  GRANT  NUMBERC*) 

ONR  N00014-75-^ 

C-0456 

'S.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

Dept,  of  Aerospace  and  Ocean  Engineering 

Virginia  Tech.  / J 

Blacksburg,  Virginia  24061 

10.  PROGRAM  ELEMENT.  PROJECT,  TASK 
AREA  A WORK  UNIT  NUMBERS 

II.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

Office  of  Naval  Research 

Fluid  Dynamics  Branch,  Code  438 

Arlington,  Virginia  Attn:  Mr.  Mort  Coooer 

12.  report  DATE 

April  1978 

IS.  NUMBER  OF  PACES 

22 

14.  monitoring  agency  NAME  S AOORESSCI/ dlfforwif  from  Coninlling  OHIe») 

IS.  SECURITY  CLASS,  (ol  Ihio  roport) 

Unclassified 

ISo.  DECLASSIFICATION/ DOWNGRADING 
SCHEDULE 

16.  DISTRIBUTION  STATEMENT  (oi  tfila  RaporO 

Distribution  unlimited 

17.  DISTRIBUTION  STATEMENT  (ol  Mia  abatraet  tnUnd  In  Block  20,  It  dlllonat  fraai  Roport) 

IS.  SUPPLEMENTARY  NOTES 

It.  KEY  WORDS  fComNnua  on  ravaraa  oldo  It  nocoooorr  and  Idontllr  hr  kloek  maabarj 

Transonic  Shocks 

Viscous-Inviscid  Interaction 

Cryogenic  Test  Facilities 

abstract  (Cpnllnum  on  rovora#  altfo  It  nocoaoorr  m4  Idttlltr  Sp  SfocA  mmbf) 

'^he  role  of  the  basic  similitude  parameters  gover 
turbulent  boundary  layer  interaction  effects  in  c 
is  studied  theoretically  for  the  non-separating  c 
Reynolds  number,  these  parameters  are  the  wall  to 
specific  heat  ratio  ^ viscosity- temperature  expo 
The  results  show  that  lack  of  temperature  ratio  s 
adverse  effect  on  interactive  skin  friction  and  h 

Cjj(v»<vyA.^"y 

ning  transonic  normal  shock- 
ryogenic  wind  tunnel  tests 
ase.  Besides  Mach  and 
total  temperatures  ratio, 
nent  and  Prandtl  number, 
imulation  has  a significantly 

@nC6  SSPardLlOn  OnScL  ^ ^ 

DO  1473  EDITION  OF  I NOV  SB  IS  OBSOLETE 


SECURITY  CLASSIFICATION  OF  THIS  RAGE  (Whmt  OMa  CnfataO 


